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O Abstract 

We study the relationship between local structural ordering and dynamical heterogeneities in a model glass-forming liquid, 
^ the Wahnstrom mixture. A novel cluster-based approach is used to detect local energy minimum polyhedral clusters and local 

a crystalline environments. A structure-specific time correlation function is then devised to determine their temporal stability. 

For our system, the lifetime correlation function for icosahedral clusters decays far slower than for those of similarly sized 
F^ but topologically distinct clusters. Upon cooling, the icosahedra form domains of increasing size and their lifetime increases 

^ with the size of the domains. Furthermore, these long-lived domains lower the mobility of neighboring particles. These 

Q structured domains show correlations with the slow regions of the dynamical heterogeneities that form on cooling towards the 

^ glass transition. Although icosahedral clusters with a particular composition and arrangement of large and small particles are 

' ' structural elements of the crystal, we find that most icosahedral clusters lack such order in composition and arrangement and 
thus local crystalline ordering makes only a limited contribution to this process. Finally, we characterize the spatial correlation 
of the domains of icosahedra by two structural correlation lengths and compare them with the four-point dynamic correlation 
length. All the length scales increase upon cooling, but in different ways. 
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I. INTRODUCTION 

The nature of the rapid increase in viscosity as hquids are cooled toward the glass transition is the subject of many 
theoretical approaches, however, there is no consensus on its fundamental mechanisnP^. One plausible scenario is that 
dense packing leads to self-induced memory effects, which causes slow dynamics^. However, the recent discovery of 
dynamic heterogeneities, i.e., spatial heterogeneities in the relaxation dynamics that emerge on supercooling^H^, is 
suggestive of the importance of a growing dynamic length scale in the slowing down approaching the glass transition. 
In addition to this dynamical phenomenon, the idea of a structural change leading to vitrification has a long histor}!^ 
and recently it has become clear that a range of glass formers show a change in structure upon the emergence of slow 
dynamics. 

Two generic types of structure have been identified: spatially extendable crystal-like ordering^Hni^ ^^d non- 
extendable polyhedral ordering^^Hi^. The latter concerns particles organized into polyhedrsPEni which cannot tile 
Euclidean space due to geometrical frustration^ ^^. Thus they form ramified structures whose fractal dimension is less 
that the dimensionality of the system, i.e., non-extendable ordering. Metallic glasses can exhibit this second kind of 
ordering^^ ^^. Other "order-agnostic" schemes have been devised where structural correlations not related to a specific 
motif have also been identified^^ ^^. 

To strengthen the link between structure and dynamics, two approaches have been employed: firstly, dynamically 
slow regions have been correlated with local ordering. Secondly, emergence of structural and dynamic length scales, 
which grow similarly, has been sought. The former links structure to dynamical heterogene ity and has been identified 
in glass formers with both crystal-like orderin^^K^^ and non-extendable polyhedral ordering^^E^'^. The latter case of 
growing structural and dynamic length scales is more controversial. One of us has identified a direct correspondence 
between the growing dynamical length scale and a structural length scale in polydisperse hard spheres with crystal-like 
ordering^EHI^ while others have not^^. In the case of binary Lennard- Jones (Kob- Andersen) mixtures, only weakly 
growing structural length scales have been identified by a "point-to-set" analysis^, while other approaches using 
static perturbation of inherent structures'^ and finite size scaling^^ find a much stronger increase of static length 
scales. Similarly, for binary hard and soft sphere mixtures, no local ordering, nor one-to-one correspondence between 
a growing structural length scale and the growing dynamical length scale has been found^^-^^-^^-^, but again static 
perturbation did suggest a growing structural length scale^'. 

In some glass for mers, local polyhedral ordering in the form of icosahedra has been related to an underlying (quasi-) 
crystalline structure^^EH. The Wahnstrom mixture^^ used here is a commonly studied binary Lennard- Jones model 
glass-former of large (A) and small {B) particles. This model can crystallize upon formation of a structure that 
contains icosahedra composed of six A and six B species around a central B-type particle, which can tile space with 
the inclusion of Frank-Kasper bonds^^. For this mo del it was shown that icosahedral order develops upon cooling and 
a link to slow dynamics has been suggestecP^'^^'^. Thus the Wahnstrom mixture provides an intriguing case where 
icosahedra are found both the in the spatially extendable crystal and in non-extendable networks. In this work, we 
find the latter. 

The purpose of this work is threefold: (i) To identify any structural origin for dynamic heterogeneities directly 
by measuring a structure-specific time correlation function of clusters detected by a topological algorithm, (ii) To 
elucidate the microscopic mechanism by which structured icosahedral domains determine the dynamics, and identify 
the role that local crystalline ordering plays in this process, (iii) To understand the behavior of correlation lengths 
for the sizes of the structured domains and those for the dynamical heterogeneities. To address these points, we use 
the Wahnstrom mixture^. 

This paper is organized as follows. In section [llj we provide details of our simulations and analysis, outlining the 
dynamic topological cluster classification used to elucidate long-lived structural motifs. In our results section |III[ 
we show that, of the structures we consider, icosahedral clusters of 13 particles are exceptionally long-lived. These 
long-lived icosahedra form a network which is correlated with dynamically slow regions, and which acts to slow down 
neighboring particles. We further show that, although the crystal structure of the Wahnstrom mixture is also formed 
of icosahedra of a specific composition of the A and B particles, these form only a small subset of icosahedra in the 
supercooled liquid. Finally we show that, of two structural correlation lengths based on these icosahedra, neither 
scales with the dynamic correlation length ^4(T). In Section IV we summarize our work. 
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Figure 1: The self- intermediate scattering function Fg (kp^t). 



II. METHODOLOGY 



Model and Simulation Details 



In the Wahnstrom mixture^, the two species of Lennard- Jones particles interact with a pair- wise potential, U{r) = 
4eQ,/3[(^^^)^^ — (^^^)^], where a and /3 denote the atom types A and B^ and Vij is the separation. The energy, 

length and mass values are caa = ^ab = ^bb^ cfbb/ctaa = 5/6, cryiB/c'"AA = 11/12 and ttia = 2771^ respectively. 
We employ NVE molecular dynamics simulations in 3D with N = 10976 particles, Na = Nb^ at constant density 
p = 1.296. Lengths, temperatures and times are quoted in units of aAA^ ^aa/^b and (mAcr^^/cAA)^^^ respectively. 
The molecular dynamics equations of motion are integrated using the Velocity Verlet algorithm with time step 



At = 0.001. The potentials are truncated at r\ 



AA 



^AB 



^BB 



2.5 and smoothed using the Stoddard- Ford method 



to ensure continuous forces^^. 

The a-relaxation time, r^, of each state point is defined using the self-intermediate scattering function (ISF) of the 

yl-type particles, F/(k, t) = wZj^i exp[ik • (rj(t) — rj(0))] V where j indexes the A-species. The angularly averaged 

F^{kp^t)^ where wavenumber kp corresponds to the first peak in the AA-partial structure factor Saa{]^)i was fitted 
with the Kohlrausch-Wilhams- Watts (KWW) stretched exponential, F^{kp,t) ^ Cexp[— (t/r^)^]. The stretching 
exponent /3 and Debye- Waller factor C are fitting parameters. Example ISFs are shown in Fig. [l] 

The highest temperature state point (T = 5.0) was initialized from a random configuration. Equilibration took 
place in the A^^T-ensemble for ^ lOOr^ using the Nose-Poincare thermostat with coupling constant Q = l.G^^, before 
further equilibration in the A^F^-ensemble for lOOOr^. On completion of the equilibration process, trajectories of 
length 300r^ were generated for analysis. The state points for T < 5.0 were obtained via a step- wise cooling process. 
Each new liquid state point was obtained by quenching instantaneously to a lower temperature from an equilibrated 
configuration of the previous higher temperature state point. An equilibration process identical to that for the T = 5.0 
liquid was performed following the quench. 

The stability of the system in the long simulation runs was checked by examining the time evolution of the ISF 
and the partial radial distribution functions ^aa(^), QABij) and gBBij). These quantities do not change throughout 
the simulation runs. There was also no drift in the numbers of the clusters detected by the TCC algorithm through 
time. Crystallization was not seen in any of the supercooled samples. 



B. The Topological Cluster Classification Method 



Groups of particles whose bond network is topologically identical to isolated cluster s whi ch minimize the potential 
energy^ are detected using the Topological Cluster Classification algorithm (TCC)P^^. While isolated clusters 
naturally experience a different environment to that of a bulk (supercooled) liquid, we note that in the case of 13 
particles, in a mean-field liquid, the icosahedron has a lower energy than local crystalline ordering^, as is the case for 
isolated clusters. The bonds are found using the modified Voronoi construction with four-membered ring parameter 
fc = 0.82. The TCC detects all 3, 4 and 5 membered shortest-path rings within the network of bonded particles, and 
from these the ground-state clusters are identified. The structure classification is independent of the composition of 
the cluster in terms of A- and B- particles. 
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Figure 2: The 13A icosahedral cluster and clusters that have icosahedral-like structure. Clusters HE and 12D occur at a- 
interlocking icosahedral sites^^, the rest are contained within a single icosahedron. Clusters in dotted rectangles are ground 
state clusters of n particles for the Wahnstrom mixture^. The colors of the particles and the bonds highlight the detection 
methods in the TCC algorithm. Grey particles are part of 3, 4 or 5 membered shortest path rings. These rings are shown by 
the white, blue, pink and green colored bonds. Yellow particles are spindle particles for the shortest path rings which together 
form the basic clusters. Red particles are additional particles bonded to a smaller clusters forming a new cluster. The detection 
routines for clusters are described in^. 
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Figure 3: Clusters with topology distinct from icosahedral order. Clusters highlighted with dotted rectangles are ground state 
clusters of n particles for the Wahnstrom mixture^^. Colors are the same as in Fig. [21 Detection routines for lOW and IIW 
clusters are described below. 
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Table I: The ground state clusters of the Wahnstrom mixture in TCC notation for n = 5 to 13 particles containing ua A-specie 
particles. The energy E is the binding energy of the cluster. Clusters highlighted in bold face are the lowest energy states for 
a given value of n. The final column is the cluster detected by the Topological Cluster Classification algorithm from the ground 
state configuration. Note that for the Wahnstrom mixture the ground state clusters with ua = and ua = ri are identical for 
each n as the AA and BB Lennard- Jones interactions are the same. The ua — clusters have therefore been omitted from the 
table for brevity. 



We define the lifetime r^ of a n-particle cluster, determined individually by the particle indices, as the time difference 
between the first and last instances of its detection by the TCC. We ensure that a cluster is detected by the TCC in 
each T^ time window during its lifetime r^ and no subset of the particles become unbonded from the others in the 
cluster during r^. The lifetime correlation function P(r£ > t) of each type of cluster is the fraction of instances of the 
cluster that have lifetime r^ greater than t. 

The clusters that are identified are all those that are the ground state clusters for the Lennard- Jones potential, 
the Wahnstrom mixture for n = 5 to 13 particles, and a range of others that are either ground states of Morse 
potentials with various interacti on le ngths, or are clusters found in face-centered cubic (FCC), and hexagonal-close 
packed (HCP) crystal structures'^ ^^. We follow the nomenclature of Doye et al.^'^ and divide the clusters into two 
categories: (i) those that are icosahedral-like in structure, i.e. clusters that are subsets of an icosahedral cluster or 
occur between sites of a- and /3-percolating icosahedral clusters^, and (ii) others with structurally distinct order. 
The clusters divided between these two categories are depicted in Fig. [2] and Fig. |3] respectively. 

For each cluster size tti, there are m-\-l different compositions for a cluster in terms of the number of A and 5-species. 
The ground-state configurations for each composition may be identified using the publicly available minimization 
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package GMIN '^ . While the clusters with the minimum binding energy E for each size n were identified in 
list ground state clusters for the other combinations of n, ua in table [ij 

The detection of structures by TCC is independent of the type of species present in the cluster, except for when 
icosahedral clusters compatible with the crystal phase are considered below. Table [l| demonstrates that the structure 
of the ground state clusters of the Wahnstrom mixture for each size n is only weakly dependent on the composition of 
the particles. The TCC algorithm presented irpSESl ^g qJ^\q ^^ identify the structure of all but two clusters in terms of 
the Morse ground state clusters. The additional two clusters lOW and IIW are specific to the Wahnstrom mixture. 

The cluster lOW is composed of six five-membered rings of particles forming a shell around a central particle. All 
particles in the shell are bonded to the central particle, and the central particle is the only nearest neighbor shared 
by all particles in each of the five-membered rings. This is the ground state for the scenario n = 10 and n^ = 9. For 
n = 11 there are multiple compositions where the IIW cluster is the ground state. This is identified as a lOB cluster 
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Figure 4: (a) Non-Arrhenius increase in structural relaxation time r^ versus the inverse of temperature. The data are fitted with 
an Arrhenius exponential at high temperatures and VFT at low temperatures. The onset temperature of the super- Arrhenius 
behavior is T* = 1.48 (vertical dashed line in both (a) and (b)). (b) The fraction of particles detected within icosahedral 13A 
clusters of all lifetimes A/'i3a(t£ > 0)/^- Note the fast increase in NisAiri > 0)/N starts around T*. 

with a single additional particle bonded to the central lOB particle. The coordination number is 10. The crystal 
structure of the Wahnstrom mixture contains icosahedral clusters and Frank-Kasper clusters of a specific composition. 
Their relationship to the icosahedral domains is discussed in the next section. 

III. RESULTS AND DISCUSSION 

A. Overall dynamics 

First we show the temperature dependence of the structural relaxation time r^ in Fig. Wa). As is usual for fragile 
glass- forming systems, we see two regimes for the behavior of the structural relaxation time delimited by an onset 
temperature for slow dynamics T* = 1.48. For T > T* an Arrhenius form, r^ = rooexp[£^oo/^B^], where Too and 
Eqo are fitting parameters, describes t^{T) well. For T < T* the Vogel-Fulcher-Tammann (VFT) form is used, 
T^ = T^ exp[DTo/{T — To)], where D is the fragility index and Tq = 0.40, for the non- Arrhenius regime of increasing 
relaxation times. Hereafter we seek structural features responsible for the super- Arrhenius behavior. 



B. Topological Cluster Classification analysis 

The temperature dependence of the number density of icosahedral structures is shown in Fig. ^h). The fraction 
of particles detected as members of icosahedral clusters is NisAi'Ti ^ 0)/A/'. Here r^ > emphasizes that here we 
consider all icosahedra rather than the time-correlated structure considered below. At high temperatures there are 
few icosahedral clusters. On cooling there is a rapid increase in Nisa{^£ ^ 0)/A^. This coincides with the onset of the 
super- Arrhenius increase in r^ [see Fig. [4[a)], suggesting a link between structure and dynamics in this fragile system. 
For the lowest temperature sample, T = 0.604, approximately 30% of particles are found to be within icosahedral 
clusters. 

In Fig. |5|a) the cluster lifetime correlation functions F{r£ > t) are shown for a range of clusters detected by the 
TCC in the lowest temperature sample (T = 0.604). The lifetime correlation function for the n = 13 icosahedral 
clusters (black line) decays much more slowly than P(r£ > t) for all the structurally distinct clusters detected by the 
TCC (grey shaded region). 

The long-time tail of the icosahedral correlation function indicates that some of these clusters preserve their local 
structure on timescales far longer than Ta- This effect may be related to previous observations that some icosahedra 
(those in which the central particle was previously in an icosahedron) are more stable than average^^. As we shall 
see below, this effect is enhanced when the icosahedra group into domains. The lifetimes of all cluster types hold no 
simple monotonic relationship to their size and frequency of occurrence. For example the four- fold symmetric n = Q 
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Figure 5: (a) Cluster lifetime correlation functions P(r£ > t) for icosahedral (denoted 13A - solid line), octahedral (6A - 
dashed line), and HCP13 (dotted line) at T = 0.604. All clusters with structure distinct from icosahedral order fall in the grey 
shaded region between P(r£ > t) for 6 A and HCP13 • (b) The fraction of particles detected within icosahedral clusters with 
lifetime t^ > t (T — 0.604). (c) The mean-square displacement (Sr'^{t)) of all particles (black line) and particles starting within 
icosahedral clusters with lifetime r^ > t as specified in the legend (colored lines - T = 0.604). 



particle octahedral cluster is smaller and occurs 8 times more frequently than icosahedral clusters, yet displays faster 
decay of P(r£ > t). Moreover, there is no trend in the rate of decay of the lifetime correlation function F{r£ > t) for the 
ground-state clusters, even though it was once conjectured that the ground-state clusters would be most stable^. For 
example, the ground state clusters 6 A, 8A and IIW have correlation functions P(r£ > t) that decay more quickly than 
the ground-state clusters that are subsets of the icosahedron (5A, 7A, SB, 9B, lOB, 12B and 13A). The correlation 
functions P(r£ > t) for the 6 A, 8 A and 11 W clusters fall within the shaded region in Fig. [sja). 

The fast initial drop in P(r^ > t) reflects the existence of large numbers of icosahedra with lifetime T£ <^ r^. The 
lifetime of these clusters is comparable to the timescale for the beta-relaxation regime where the particles fluctuate 
within their 'cage' of neighbors. These clusters are not representative of the underlying viscous structure of the liquid 



that exists on the timescales of the dynamic heterogeneities ^ r^. Fig. |5[b) demonstrates that particles within these 
short-lived icosahedral clusters make only a small additional contribution to the total fraction of particles detected 
within all icosahedral clusters. Icosahedral clusters overlap and particles can be members of multiple icosahedral 
clusters simultaneously, and the results of Fig. [5|b) indicate that short- and long-lived icosahedra mainly lie in the 
same regions of the liquid. Here we do not distinguish the short- and long-lived icosahedral clusters structurally. 

The mean-square displacement (MSD), {Sr'^{t)) — l/A/"^^ ki(^) ~ ^i^y? is shown in Fig. p[c) for all particles 
(solid black line) and for particles initially within icosahedral clusters (colored lines). The dynamics of particles 
within icosahedra of all lifetimes are slower than the system- wide average (c.f. black line). The longer the lifetime of 
the icosahedra, the slower the particles become. Since some icosahedra last for very long terms [Fig. [5|b)], we expect 
that these particles may exhibit very low mobilities. 



C. Network of icosahedra 



The domains of icosahedra that form on cooling are shown in Fig. [6J At high temperature icosahedra are predom- 
inantly isolated [Fig. [6]^ a)]. Upon cooling, the icosahedra overlap and join together [Fig. |6jb)] to form (transient) 
networks at low temperatures [Fig. [6]^c)]. In order to investigate the structure of the domains, we calculate the radius 
of gyration Rg^ R^ = ^^ ^^ (r^ — r^)^, where m is the number of particles in the domain and the sum extends 
over all pairs of particles in the domain. The radius of gyration shows power-law growth with an exponent of 0.47 
with the size of the domain m [Fig. [oFd)]. In other words the domains have a fractal dimension df c^ 2, indicating 
that they are non-space-filling. We find that the individual icosahedra have enhanced stability as the size of the 
icosahedral domains grow [Fig. [oFe)]. The time r^ is the mean lifetime of icosahedra taken from size m domains in 
each configuration. For T = 0.620, this grows almost five-fold as the domain size increases from isolated icosahedra 
(m = 13) to extended clusters of icosahedra (m ^ 200). For higher m, there is some suggestion that ri may saturate. 
At present, it is hard to be sure. Our data are limited by system size, so much larger simulations would be required 
to properly address this issue. If indeed, there is saturation in r^ with 772, this would suggest a maximal stability of 
clusters. Thus that upon percolation where m diverges, f£ need no increase further and full dynamical arrest would 
not be expected, as we indeed find. 

We now consider the effect this network has on the remainder of the system. In Fig. |6Ff), we plot the mean-square 
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Figure 6: (a)-(c) Domains of icosahedral clusters form on cooling from high to low temperature (slices through 3D simulation 
box). Particles in icosahedral clusters are shown full size in green, other particles are blue dots, (d) The radius of gyration Rg 
of the domains versus the number of particles in the domain ra for T — 0.620. Rg is well fitted by mP'^^ indicating the domains 
have a fractal dimension df ^ 2. (e) The mean lifetime of icosahedral clusters ri versus the domain size m. (f) Icosahedral 
domains retard the motion of neighboring particles. The MSD {6r'^{Th)) as a function of distance from icosahedral domains d 
(solid line). The dotted line is the MSD over th of all particles not in icosahedra, independent of d. Note that the particles 
located around d = 1 amount to 40% of the system at this temperature. 

displacement of the non-icosahedral particles (^r^(r^)) against the distance d from a domain of icosahedra at time 
t = 0. Here the timescale is taken as the maximum of the four-point dynamic susceptibility X4(^W ^h — ^a- The first 
nearest neighbors of the domains have suppressed mobility, due to coupling to the stable and slow domains. The second 
and third shells of neighbors of the domains have higher mobility. This connection between the domains of icosahedra 
and the spatially heterogeneous dynamics is depicted in Fig. [TFa). The transient network of slow icosahedral particles 
retards the dynamics of neighboring particles, leading to a mechanism for the formation of dynamical heterogeneities. 
This is reminiscent of a point-to-set mechanism^^-^ whereby pinned (immobilized) particles inhibit the motion of 
their surroundings. Of course here the particles in long-lived icosahedra are slow, rather than immobilized, so the 
network is transient, but its neighbors are nevertheless retarded. 



D. Relationship of icosahedral domains to crystalline order 



Pedersen et al. identified that the crystal structure of the Wahnstrom mixture has the structure of the MgZn2 Laves 
phase, which is formed of 13A icosahedra and Frank-Kasper bonds with specific compositions arrangements of the A- 
and ^-species^^. They found that icosahedra within the crystal structure are composed of six of the larger A-particles 
and seven of the smaller 5-species. Furthermore of all the possible arrangements of the A- and 5-particles within 
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Figure 7: (a) The relationship between icosahedra and spatial heterogeneities in the dynamics at T = 0.604 (2D slice of 
simulation box). Mobility ((^r^(T^))) is depicted by the size of the particle. Large particles have (^r^(r^)) < 0.043, small 
particles otherwise. Particles always within icosahedral clusters with lifetime r^ > r^ are colored green, and others are grey, 
(b) Growing structural ^Rg (squares) and ^si3A (crosses), and dynamical ^4 (circles) correlation lengths. 
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Figure 8: (a) Relationship between icosahedral 13A domains and the crystalline ' Frank- K asp er' and icosahedral clusters for 
T = 0.604. The quantity Nc{ri > t)/N is the fraction of particles detected within clusters lifetime ri > t. The dashed line 
is for Frank-Kasper clusters, and the dotted line is for icosahedral clusters with a B-particles in the center and arrangements 
of A- and 5-species in the shell of the icosahedron compatible with the bulk crystal phase^^. Both of these clusters exist in 
the MgZn2 structure of the Wahnstrom crystal. The solid line is for all icosahedra, irrespective of the arrangement of A- and 
5-species. (b) Snapshot at T = 0.604 showing the overlap between icosahedral domains and Frank-Kasper clusters. Green 
particles are the icosahedral domains, red particles are those in both Frank-Kasper and icosahedral clusters, and black are 
Frank-Kasper particles not in icosahedral domains. The small blue particles are not members of icosahedral or Frank-Kasper 
clusters. 



the icosahedral cluster for this composition, only two are present in the crystal structure. They also suggested that 
Frank-Kasper clusters, consisting of two A-specie particles bonded to six common 5-specie neighbors, that occur in 
the crystalline lattice, act to stabilize the supercooled liquid. Here we discuss the relationship between the icosahedral 
domains we identify and crystalline structure in the supercooled liquid. 

In Fig. IsFa) we present data for the fraction of particles detected within clusters of lifetime ri > t for three 
types of clusters. First are for the icosahedral clusters (13 A) that form extend domains through the system. These 
clusters can be made of any number A- and 5-species, i.e. icosahedra with any compositions of A- and 5-particles 
and any arrangement of the particles within the cluster for a given composition. Second are the icosahedral clusters 
compatible with the crystal structure^. There are two different arrangements that consist of six A-particles and six 
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Figure 9: (a) The four-point dynamic susceptibility X4.{t)- The maximum of X4.{t) occurs at time r^. (b) The four-point 
dynamic structure factor SA(k^Th)- The solid lines are fits to the data with the Ornstein-Zernike equation E]. (c) The structure 
factor of the icosahedral particles Si2,A(k). The solid lines are fits to the data with the OZ function. 



5-particles arranged around a central 5-particle. These are identified by the number of connected 5-species in two 
sets on the shell of the icosahedral cluster which are separated by a ring of six A-particles. One arrangement has two 
sets of three connected 5-particles, and the other one pair and one quadruplet of 5-particles. The particles within 
crystal-compatible icosahedral clusters are necessarily a subset of the particles within all icosahedral The third cluster 
is the Frank-Kasper cluster, which consists of a bonded pair of A-species surrounded by a ring of 6 5-particles. 

For all cluster lifetimes there are more particles in icosahedral clusters than in Frank-Kasper or crystal compatible 
icosahedra. It is strictly necessary that particles in crystal compatible icosahedral clusters are fewer than in icosahedral 
clusters, as the former are by definition icosahedra themselves. However we emphasize that very few particles (< 5%) 
are identified in the clusters compatible with the crystal for the lowest temperature state point we studied. We find 
for T = 0.604 that 83% of particles in Frank-Kasper clusters also lie within the icosahedral domains. The overlap of 
particles in Frank-Kasper clusters and icosahedral domains is shown in Fig. IsFb). 

We conclude that crystalline structuring in supercooled Wahnstrom liquids is a relatively small component of 
wider icosahedral ordering. The icosahedral ordering is frequently incompatible with the crystal structure due to the 
composition, i.e. ua is frequently not equal to 6, or because the arrangements of the A- and 5-species within the 
cluster are inconsistent with the crystalline icosahedra. 



E. Comparison of structural and dynamic length scales 

Finally we investigate how the four-point correlation length for the dynamic heterogeneities ^4 compares to two 
static correlation lengths for the sizes of the domains of icosahedra. The first of these is the radius of gyration of the 
icosahedral clusters ^Rg = i?J^^((m)/13)^-^^, where (m) is the ensemble average of the domain size and Rl^^ = 0.95 
is the radius of gyration of a single icosahedral cluster. This demonstrates the growth in the domains on cooling. The 
second, which we term ^sisA which is calculated from fitting the Ornstein-Zernike (OZ) equation [Eq. El to the low 
k behavior of the static structure factor restricted to the icosahedral particles SisAJk). 

We follow Lacevic et al. in our calculation of the dynamic correlation length ^J^. We define an overlap function, 
w{\yj{{)) — r^(t)l), to be unity if |rj(0) — vi(t)\ < a, otherwise, where a = 0.3. It defines the particle sites at time 
that are occupied by a particles at time t. The fraction of overlapping particles between two configurations is 
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The four-point dynamic susceptibility Xi{t) is then given by 
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X4{t) 



N^keT 
The time at which Xi{t) has a peak is r^ [Fig. |9]ja)]. 



[{Qm - {Qim. 



(1) 



(2) 



10 



The four-point dynamic structure factor ^4 (k, t) is defined as 

54(k,t) = ^(^exp[ik-r,.(0)]w;(|r,(0)-r,(t)|)^exp[ik-r„(0)]«;(|r^(0)-r„(t)|)), (3) 

jl mn 

where j, /, tti, n are particle indices. For time r/^, the angularly averaged version is S4^{k^Th). The dynamic correlation 
length ^4 is calculated by fitting the OZ equation, 

to S4^{k^Th) for /c < 3^ [Fig. M^b)]. To calculate ^513^7 we follow a similar approach as for ^4. The structure factor 
of the icosahedral particles is given by 

where j, / index particles that are detected within an icosahedral cluster and Nisa = ^isAi'Ti ^ 0) is the total number 
of icosahedral particles in the configuration. The correlation length ^513^ is found by fitting the Ornstein-Zernike 
(OZ) equation to the low k values of the radially averaged icosahedral structure factor SisA{k) [Fig. |9|c)]. 

The temperature behavior of these lengths is shown in Fig. W[h). All three correlation lengths increase on cooling. 
For the lowest temperature studied, the network of icosahedra percolates, thus ^ng cannot be defined, suggesting a 
divergence in a structural length scale at a temperature higher than that inferred from the dynamic correlation length 
^4. Thus percolating domains of icosahedra do not imply structural arrest (see the discussion in Sec. III.C): since 
icosahedra have a finite lifetime, the formation of a network does not lead to arrest, unlike the case of colloidal gels^. 
The length ^5'i3a grows more slowly than ^4 as the domains of icosahedra are smaller than the immobile regions 
selected by the four-point dynamic susceptibility X4(^) [see Fig. ufa)]. This result is expected because domains of 
icosahedra in this liquid form rarefied networks with df c^ 2 which do not fill space. It is unlikely that this difference 
would be rectified by including the immediate neighbors of icosahedral domains when calculating the structural 
correlation lengths, because there are large numbers of neighbors of icosahedral domains that are found at the lowest 
temperatures. We also find that there are a fraction of the particles within icosahedral clusters that are not selected in 
the calculation of ^4 (i.e. their displacement over Th is greater than 0.3, meaning that there is not a perfect correlation 
between the structured and the slow particles. Note X4 has been shown to exhibit dependence upon system size for 
A^ ^ lOOG^^. We expect that such effects are reasonably small here, and have taken care to only consider temperatures 
where all our measured lengths are much smaller than the system size. While system size effects cannot be ruled out, 
we do not believe these make a significant impact on the conclusions we draw. 

The result for this Wahnstrom mixture contrasts with liquids with crystal-like order. In some work, this has been 
found to scale with ^4^^^^^. Other methods have also shown that structural correlation lengths are rather larger in 
the case of such liquids^. We emphasize that the icosahedral clusters which form the domains have arrangements 
of the A- and 5-species that are rarely compatible with the crystal structure. Only 10% of the particles within 
icosahedral domains are members of clusters in the crystal structure. Thus, although an icosahedral structure itself is 
a key structural motif of the Frank-Kasper crystal, icosahedral clusters formed in a supercooled liquid have a rather 
weak link to the crystal structure [Fig. |8[b)], presumably due to entropic frustration-- of the chemical arrangement 
of the species within the cluster. 

The difference in the nature of the structural fiuctuations between this model and those with spatially extendable 
crystal-like order may be responsible for the difference in the T- depend ence of ^4: ^4 ex ( T — To )~^-^^ is found here 
whereas ^4 ex (T — Tq)"^/^ for liquids with extendable crystal-like^SHni or amorphous ordei^^^^. Thus, the question 
remains as to the most appropriate structural correlation length to explain the viscous slowing down of supercooled 
liquids, and how order- speci fic correlation lengths, such as ^Rg and £,si3A^ might relate to order-agnostic structural 
correlation lengths^^^^^^^^ESl. We also note that Fig. p[f) strongly indicates that the effect of the domains on the 
surrounding particles extends around a from the domains, showing that the dynamical effect of icosahedral clusters 
is not just limited to their domains. Thus icosahedral clusters may act to slow some of their neighboring particles to 
form regions of slow dynamics, which appear correlated with spatial patterns of dynamic heterogeneities [Fig. 4(a)], 
although this should be carefully checked. 
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IV. SUMMARY 

We have demonstrated that studying the hfetimes of structural ordering within the Wahnstrom supercooled liquid 
can be used to directly identify the order that constitutes the slow regions of dynamic heterogeneities. The structures 
responsible, the icosahedra, combine to form non-space filling domains on cooling that become more stable as their 
size increases. Here, the crystal structure also has an icosahedral component, but icosahedra compatible with the 
crystal are rarely found in the supercooled liquid. Particles that neighbor the domains have suppressed mobility. 
These results imply that growing domains of icosahedra play an important in the microscopic origins of dynamic 
heterogeneities in this system. At the same time, the lack of a direct link in the length scales associated with 
structure and dynamics requires further study on the origin of slow dynamics including the mechanism of propagation 
of slowness from icosahedral clusters to their neighbors. 

AcknoAvledgements 

We thank M. Leocmach for his kind help in bond orientational order analysis and its comparison with TCC 
analysis. A.M. is funded by EPSRC grant code EP/E501214/1. C.P.R. thanks the Royal Society for funding. H.T. 
acknowledges support from a grant-in-aid from the Ministry of Education, Culture, Sports, Science and Technology, 
Japan and the Aihara Project, the FIRST program from JSPS, initiated by CSTP. This work was carried out using 
the computational facilities of the Advanced Computing Research Centre, University of Bristol. 



* Electronic address: p addy.royall@bristol.ac.uk| 

^ Electronic address: tanaka@iis.u-tokyo.ac.jp 

^ L. Berthier and G. Biroli, Rev. Mod. Phys. 83, 587 (2011). 

^ W. Gotze, Complex Dynamics of Glass-Forming Liquids: A Mode-Coupling Theory (Oxford University Press, Oxford, 2008). 

^ M. M. Hurley and P. Harrowell, Phys. Rev. E. 52, 1694 (1995). 

^ M. Ediger, Annu. Rev. Phys. Chem. 51, 99 (2000). 

^ L. Berthier, G. Biroli, J. P. Bouchaud, L. Cipelletti, and W. Van Saarloos, Dynamical heterogeneities in glasses, colloids, 
and granular media (Oxford University Press, Oxford, 2011). 

^ F. C. Frank, Proc. Roy Soc. A 215, 43 (1952). 

^ H. Shintani and H. Tanaka, Nat. Phys. 2, 200 (2006). 

^ H. Tanaka, T. Kawasaki, H. Shintani, and K. Watanabe, Nat. Mater. 9, 324 (2010). 

^ F. Sausset and G. Tarjus, Physical Review Letters 104, 065701 (2010). 
^° M. Leocmach and H. Tanaka, Nat. Comm. 3, 974 (2012). 
^^ H. Jonsson and H. Andersen, Phys. Rev. Lett. 60, 2295 (1988). 
^^ T. Kondo and K. Tsumuraya, J. Chem. Phys. 94, 8220 (1991). 
^^ T. Tomida and T. Egami, Phys. Rev. B 52, 3290 (1995). 
^^ R. Jullien, P. Jund, and P. Caprion, Phys. Rev. E 54, 6035 (1996). 

^^ M. Dzugutov, S. L Simdyankin, and F. H. M. Zetterling, Phys. Rev. Lett. 89, 195701 (2002). 
^^ E. Lerner, L Procaccia, and J. Zylberg, Phys. Rev. Lett. 102, 125701 (2010). 
^^ U. R. Pedersen, T. B. Schr0der, J. C. Dyre, and P. Harrowell, Phys. Rev. Lett. 104, 1 (2010). 
^^ D. Coslovich, Phys. Rev. E 83, 8 (2011). 

^^ G. Tarjus, S. A. Kivelson, Z. Nussinov, and P. Viot, J. Phys.: Condens. Matter 17, R1143 (2005). 
^° T. Schenk, D. Holland-Moritz, V. Simonet, R. Bellissent, and D. M. Herlach, Phys. Rev. Lett. 89, 75507 (2002). 
2^ D. B. Miracle, Nature Materials 3, 697 (2004). 

22 G. Biroli, J. P. Bouchaud, A. Cavagna, T. S. Grigera, and P. Verrochio, Nat. Phys. 4, 771 (2008). 

23 M. Mosayebi, E. Del Gado, P. Ilg, and H. C. Ottinger, Phys. Rev. Lett. 104, 205704 (2010). 
2^ F. Sausset and D. Levine, Phys. Rev. Lett. 107, 045501 (2011). 

2^ W. Kob, S. Roldan- Vargas, and L. Berthier, Nat. Phys. 8, 164 (2011). 

2^ C. Cammarota and G. BiroU, EuroPhys. Lett. 98, 36005 (2012). 

2'^ G. M. Hocky, T. E. Markland, and D. R. Reichman, Phys. Rev. Lett. 108, 225506 (2012). 

28 A. J. Dunleavy, K. Wiesner, and C. P. Royall, Phys. Rev. E 86, 041505 (2012). 

29 B. Charbonneau, P. Charbonneau, and G. Tarjus, Phys. Rev. Lett. 108, 35701 (2012). 
^° S. Karmakar and L Procaccia, arXiv: 1204.6634 (2012). 

^^ A. Widmer-Cooper and P. Harrowell, Phys. Rev. Lett. 96, 185701 (2006). 

^2 M. Mosayebi, E. Del Gado, P. Ilg, and H. C. Ottinger, J. Chem. Phys. 137, 024504 (2012). 

^^ M. Dzugutov, Phys. Rev. Lett. 70, 2924 (1993). 

^^ G. Wahnstrom, Phys. Rev. A 44, 3752 (1991). 

12 



35 
36 
37 

38 
39 

40 
41 
42 
43 
44 

45 
46 
47 
48 



S. Stoddard and J. Ford, Phys. Rev. A 8, 1504 (1973). 

S. Nose, J. Phys. Soc. Jpn. 70, 75 (2001). 

J. Doye and L. Meyer, Phys. Rev. Lett. 95, 1 (2005), ISSN 0031-9007. 

S. R. Wmiams, arXiv:0705.0203 (2007). 

C. P. Royall, S. R. Wilhams, T. Ohtsuka, and H. Tanaka, Nat. Mater. 7, 556 (2008). 
S. Mossa and G. Tarjus, J. Chem. Phys. 119, 8069 (2003). 

A. Mahns, Ph.D. thesis. University of Bristol (2013). 

J. P. K. Doye, D. J. Wales, and R. S. Berry, J. Chem. Phys. 103, 4234 (1995). 

D. J. Wales and J. P. K. Doye, J. Phys. Chem. A 101, 5111 (1997). 

Gmin: A program for finding global minima and calculating thermodynamic properties from basin- sampling, URL http: 

//www- wales . ch . CcLin . ac . uk/GMIN/ 

The Cambridge cluster database, URL http://www-wales.ch.cain.ac.uk/CCD.htinll 

C. Cammarota and G. Biroli, Proc. Natl. Acad. SciTU.S.A. 1097^850 (20T2)^ 

N. Lacevic, F. W. Starr, T. B. Schr0der, and S. C. Glotzer, J. Chem. Phys. 119, 1605094 (2003). 

S. Karmakar, C. Dasgupta, and S. Sastry, Proc. Nat. Acad. Sci. U.S.A. 106, 3675 (2009). 



13 



